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Abstract : 

A  cylindrical  water  model  was  used  to  visualize  the  flow  field  in 
a  jet-stirred  bath.  The  model  was  made  of  transparent  plexiglass.  A 
plane  of  light  illuminated  a  cross  section  of  the  cylinder  containing 
its  axis  and  the  point  of  injection;  and  thereby  illuminating  the  small 
alumina  particles  suspended  in  the  water.  The  flow  field  was  photo¬ 
graphed,  and  dye  visualization  studies  were  also  conducted. 

A  simple  modification  to  the  TKACH  computer  code  was  made  so  as 
to  account  for  the  geometry  of  the  model.  Computer  runs  were  made  for 
different  injector  Reynolds  numbers. 

Conclusions  were  made  regarding  the  mean  and  fluctuating  nature 
of  the  flow  field. 
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CHAPTER  I 


INTRODUCTION 


1.1  Statement  of  the  Problem 

This  study  is  aimed  at  obtaining  an  understanding  of  the  recir¬ 
culating  flow  field  in  a  modified  single  injector  reactor  (SIR).  This 
consisted  of  an  injector  stirred  liquid  hath  in  a  cylinder  modified  to 
have  outflow  on  both  ends.  The  primary  objective  is  to  qualitatively 
obtain  a  picture  of  the  flow  structure  as  a  function  of  the  injector 
Remolds  number  and  the  ratio  of  outflow  at  the  two  ends  of  the  model. 


1.2  Origin  and  Relevance  of  the  Sttidv 

The  flow  field  and  heat  transfer  characteristics  of  a  iet  stirred 
bath  arc  of  significant  interest.  Such  a  jet  mixing  process  is  employed 
in  many  combustion  applications.  An  example  of  such  a  system  is  a 
boiler  reactor.  In  this  case,  the  fuel  is  a  molten  Lithium  bath  and 
the  oxidizer  is  injected  into  the  reactor.  The  reaction  is  exothermic, 
and  the  resulting  heat  generated  is  transferred  to  a  boiler  which  provides 
steam  to  the  power  cycle,  ihc  jet  induced  mixing  causes  an  augmented 
heat  transfer  rate  to  the  water. 

Work  has  been  done  to  obtain  empirical  relationships  for  the  heat 
transfer  coefficients  (1),  and  to  develop  analytic  models  for  the 
velocity  distribution  and  heat  transfer  coefficients.  However,  it  was 
felt  that  in  order  to  develop  suitable  predictor  techniques  for  such 
situations,  a  more  fundamental  understanding  of  the  flow  structure  is 
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1.8  I' rev  ions  Related  Studios 

Mixing  in  ducted  turbulent  jets  was  examined  by  Becker,  llottel 
and  Williams  (2).  They  studied  the  effect  of  the  Craya-Curtet 
similarity  parameter  on  turbulent  mixing  patterns  in  axisvmmetric , 
constant  density  confined  iet  flows,  in  the  regime  in  which  recircu¬ 
lation  occurred.  Their  experimental  set-up  consisted  of  a  round  iet 
discharging  axially  into  a  cylindrical  duct  fed  by  a  uniform  stream 
who^e  entrance  momentum  relative  to  that  of  the  iet  varied  from  zero 
to  moderately  high  values.  The  fields  of  mean  velocity  and  mean  con¬ 
centration  were  mapped  and  analyzed.  They  concluded  that  at  small 
values  of  the  Crava-Curtct  number  (Ct),  the  initial  stream  flow  into 
the  duct  falls  short  of  the  entrainment  needed  by  the  iet.  The 
deficiency  is  made  up  by  fluid  recirculated  from  downstream.  They 
describe  properties  of  the  recirculation  eddy  and  the  post  iet  mixing 
zone  downstream  of  reattachment.  It  should  be  noted  that  in  the  SIR 
the  initial  stream  velocity  is  zero  and  (it  approaches  cero  implying 
total  recirculation. 

Barchilon  and  Curtet  (S')  in  attempting  to  improve  the  assumptions 
made  in  the  approximate  theory  of  confined  jets  (4,S,6)  made  some 
interesting  visualization  studies.  They  detected  unusually  high  turbu¬ 
lence  levels  in  the  iet  and  recirculation  regions.  They  found  that  the 
recirculation  eddy  which  initially  appeared  at  a  Craya-Curtet  number  of 
O.P7('  extended  itself  as  Ct  diminished  and  at  Ct  =  0.075  it  occupied  . 


the  complete  duct.  Tow  Craya-Curtet  numbers  are  indicative  of  high 
turbulence  and  large  eddies  implying,  better  mixing.  There  have  been  many 
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other  studies  on  recirculating  flows  in  ducted  jets,  they  have  been 
well  summarized  in  (1),  (2)  and  (3). 

Abramovich  (7)  developed  a  theory  for  predicting  flow  patterns  in 
a  dead  ended  channel  stirred  by  a  submerged  injector,  which  is  the 
situation  existing  in  the  SIR.  Abramovich  treated  the  flow  region  as 
two  separate  parts;  the  first  in  which  the  turbulent  jet  spreads  through 
a  counterflowing  stream  of  fluid,  and  the  second  in  which  the  streamlines 
turn  in  accordance  with  the  laws  of  motion  of  an  ideal  fluid.  Solu¬ 
tions  for  the  velocity  fields  were  obtained  by  solving  the  conservation 
equations  in  each  area.  Only  incompressible  fluids  were  considered 
and  corrections  were  anplied  to  these  solutions  to  account  for  non- 
uniform  exit  velocity  profile.  Comparison  with  experimental  results 
showed  good  agreement  of  the  predict i on  in  the  first  region  but  very 
poor  prediction  in  the  second  region. 

Some  preliminary  work  at  AR1.  analytically  determined  the  velocity 
distribution  and  the  wall  heat  transfer  coefficient  for  a  closed  ended 
cylinder  similar  to  the  SIR.  The  model  divides  the  flow  field  into 
four  parts:  (i)  Jet  with  core  region,  (ii)  Developing  or  spreading 
jet  reg. ion,  (iii)  Viscous  dissipation  region,  and  (iv)  Stagnation  region 
The  analysis  employs  previous  free  jet  (f?)  and  open  ended  confined  jets 
(?,!')  to  develop  equations  which  are  solved  numerically.  A  closed  end 
effect  is  generated  by  incorporating,  a  stagnation  region.  The  fluid 
was  assumed  to  he  incompressible  and  both  local  and  overall  heat  trans^c 
calculations  made.  No  comparison  to  experimental  data  was  made. 
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Kemey,  ct  al.  (10)  and  Keimer,  et  al.  (11)  studied  the  penetra¬ 
tion  characteristics  of  vapor  jets  submerged  in  subcooled  baths. 
Experimental  results  were  used  to  develop  empirical  relationships  for 
the  penetration  length  (10).  h’eimer,  et  al.,  using  a  variable  density 
single  fluid  model  for  the  two  phase  flow,  with  the  turbulent  mixing 
process  treated  by  an  entrainment  law  found  good  correlation  with 
results  of  their  experiments  and  those  of  earlier  investigators  over 
a  wide  range  of  operating  conditions  and  injector  geometries. 

Averv  (12)  considered  the  case  of  a  gaseous  oxidizer  jet  dis¬ 
charging  through  an  injector  submerged  in  a  bath  containing  a  liquid 
metal.  The  two  phase  turbulent  combustion  process  was  analyzed  in  a 
manner  similar  to  conventional  diffusion  flames.  A  variable  density 
single  fluid  model  was  used  to  represent  the  two  phase  mixture  in 
conjunct  ion  with  an  entrainment  assumption  appropriate  for  systems  with 
large  density  differences  between  the  jet  and  ambient  fluid.  The 
theory  provided  a  general  relationship  for  estimating  combustion  and 
condensation  lengths  for  turbulent,  unconfincd,  forced  jets  based  on 
injector  and  ambient  properties.  Results  indicate  that  mixing  charac¬ 
teristics  of  gas-liquid  icts  when  corrected  for  density  variations  are 
similar  to  single  phase  _iets. 

Thomson  (1)  studied  the  heat  transfer  for  a  jet  induced  mixing 
process  in  a  cylindrical  chamber.  A  heated  fluid  was  injected  centrally 
from  one  end  of  the  chamber  into  a  liquid  bath.  Both  ends  of  the 
chamber  were  insulated  so  only  heat  transfer  through  the  side  wall 
was  considered.  Data  was  obtained  for  water  and  ethylene  glycol. 
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varying  the  injector  diameter,  fluid  temperature,  jet  momentum  and 
chamber  dimensions.  An  analysis  assuming  total  iet  momentum  dissipa¬ 
tion  into  shear  stress  along  the  chamber  wall  was  developed  using  the 
Reynold  analogy  to  determine  the  characteristic  parameters  of  the 
system,  flood  agreement  between  the  analysis  and  data  set  was  obtained. 
From  the  analysis,  the  data  was  correlated  in  terms  of  a  Stanton 
number,  Reynolds  number,  Prandtl  number  and  chamber  length  to  diameter 
rat io. 

f.osman,  et  al.  (13)  present  numerical  solutions  for  seven  dif¬ 
ferent  two  dimensional  turbulent  elliptic  flows.  The  solution 
procedure  employed  is  that  embodied  in  the  TliACll  computer  program  (Ml. 
The  calculated  properties  of  the  seven  flows  were  connnred  with 
experimental  results  which  validated  the  numerical  procedure  and  the 
two  equation  turbulence  model  described  by  l.aunder  ct  al.  (151.  Of 
particular  interest  here  is  their  investigation  of  flow  downstream 
of  an  ax i symmet r ic ,  sudden  enlargement  in  channel  flow.  The  flow- 
structure  is  only  slightly  different  from  that  in  the  SIR.  They  observe 
strong  recirculation  and  their  results  compare  well  with  the  measure¬ 
ments  of  Back  and  Roschke  (1(0  . 

There  has  also  been  a  substantial  amount  of  work  on  flow  visualiza¬ 
tion  in  mixing  vessels.  Reynold  (1")  used  a  dye  to  visualize  a  iet 
directed  into  a  water  bath,  .let  breakdown  was  observed  as  a  function 
of  Reynolds  number.  Reynold  observed  five  categories  of  flows  which 
he  classified  as  (going  from  lower  Reynold  numbers  to  higher  ones), 
shearing  miffs,  symmetric  condensation,  sinuous  undulations,  pedal 
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breakdown,  and  confused  breakdown  (at  Reynold  numbers  higher  than  300). 

McNaughton  and  Sinclair  (IS)  studied  1 iquid-into-liquid  jets  in 
short  cylindrical  vessels  using  aoueous  blue  tracer  solution  in  con¬ 
junction  with  transparent  cylinder  walls.  Four  main  types  of  jets 
were  observed:  dissipated  laminar  jets  (Reynold  numbers  less  than  300) , 
fully  laminar  jets  (Reynold  numbers  between  300  and  1000),  semi-turbulent 
(Reynold  numbers  between  1000  and  3000)  and  fully  turbulent  (Reynold 
number  greater  than  3000).  The  laminar  length  of  sub-turbulent  jets 
was  investigated  and  correlated  with  the  Reynolds  number  and  geometric 
parameters.  Crow  and  Champagne  (10)  also  used  dyes  to  study  the 
evolution  of  jet  instability  with  advancing  Reynolds  number.  They 
also  employed  schlieren  and  light  scattering  visualization  techniques 
to  study  the  jet  structure. 

Back  and  Roschke  (16)  investigated  the  flow  of  water  through  an 
abrupt  circular  channel  expansion.  The  shear  layer  between  the  central 
jet  and  reverse  flow  region  along  the  wall  downstream  behaved  differently 
in  the  observed  regimes.  Dye  injection  into  the  shear  layer  provided 
a  means  of  observing  the  jet  reattachmcnt:  lengths  in  the  various  regimes 
from  laminar  to  fully  turbulent.  The  reattachmcnt  lengths  were  found 
to  increase  from  laminar  flow  to  reach  a  maximum  and  then  decrease  for 
highly  turbulent  flows. 

White,  ct  al.  (20)  used  a  particle  visualization  technique  to 
study  the  recirculating  flow  of  a  submerged  jet  discharging  into  a 
dead  ended  vessel.  They  used  this  to  studv  the  tangential  vorticity, 
and  the  visualization  study  was  not  extensive.  A  similar  study  was 


carried  out  by  Fkhcian  and  Iloul t  (21)  to  study  the  intake  process  in 
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an  internal  combustion  engine. 

Numerous  other  authors  have  done  work  similar  to  what  has  been 
discussed  in  this  section.  However  no  paper  directs  itself  at 
obtaining  an  understanding  of  the  overall  flow  structure  in  jet  stirred 
baths . 


1.4  Scope  and  Objective  of  this  Study 

From  this  discussion,  it  is  evident  that  a  better  understanding 
of  the  flow  structure  in  the  SIR  is  needed.  The  purpose  of  this 
investigation  was  to  visualize  the  flow  structure  in  the  SIR  in  order 
to  obtain  such  qualitative  information. 

In  order  to  simulate  the  SIR,  water  was  injected  axially  into  a 
cylinder  which  provided  for  a  variable  outflow  at  both  ends  at  a  low- 
momentum  flux.  The  flow  was  seeded  with  fine  alumina  particles,  a 
high  intensity  light  source  with  a  condensing  lens  focused  a  thin 
plane  of  light  onto  a  cylinder  axis.  A  camera  located  perpendicularly 
above  this  illuminated  plane  recorded  the  flow  pattern  at  various 
inlet  momentum*  and  outflows  (at  the  two  ends).  Dye  visualization 
techniques  were  used  both  in  the  jet  and  in  the  shear  layer,  but  with 
limited  success. 

An  analysis  was  carried  out  using  a  modified  form  of  the  TF.AC1I 
comnutcr  program  (IS,  14).  Velocity,  pressure,  turbulent  kinetic 
energy  and  dissipation  profiles  were  obtained  for  the  SIR  geometry  as  a 
function  of  Reynold  number.  A  comparison  is  made  between  experimental 


and  analvtic  trends. 
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CHAPTER  II 
THEORY  AND  ANALYSIS 


2.1  Analog  Technique 

A  water  model  was  employed  to  simulate  the  flow  field  in  a  SIR, 
so  as  to  facilitate  flow  visualization.  In  order  to  conduct  such  an 
investigation  with  a  water  model,  it  is  necessary  to  satisfy  the 
similarity  conditions  (21). 

Similarity  can  be  classified  as  follows: 

1.  Geometric  similarity 

2.  Kinematic  similarity 

3.  Dynamic  similarity 

4.  Thermal  similarity 

Geometric  simi  lari  ty  requires  a  single  scale  ratio  between  critical 
dimensions  of  the  prototype  and  model.  In  the  case  of  the  SIR-model 
this  scale  ratio  is  taken  as  unity,  hence  satisfying  geometric 
simi laritv. 

Kinematic  similarity  requires  that  the  slopes  of  the  streamlines 
are  the  same  in  the  model  and  the  prototype,  implying  that  the  cor¬ 
responding  flows  are  alike.  The  assumption  that  the  fluid  in  both 
cases  are  incompressible  (assuming  the  compressibility  effects  in  the 
gaseous  oxidizer  are  negligible)  leads  to  the  conclusion  that  the 
continuity  equation  and  hence  the  streamlines  are  identical  for  the 
prototype  and  model  (22). 
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Dynamic  similarity  requires  that  the  forces  acting  on  the  f’uid 
elements  have  the  same  ratio  to  each  other  at  any  instant.  In  the 
SIR  hath  both  inertia  and  viscous  forces  are  of  importance,  and  it 
would  seem  logical  to  conclude  that  the  Reynolds  numbers  should  be 
equal  for  the  model  and  prototype.  This  would  involve  the  problem 
of  determining  the  Reynolds  number  for  a  two-phase  system.  This 
problem  is  avoided  by  considering  that  the  inertia  forces  in  the  jet 
are  of  greater  significance  (high  Reynolds  numbers)  than  the  viscous 
forces  in  the  iet.  The  injector  momentum  flux  for  the  model  was  taken 
as  the  same  momentum  flux  that  is  present  in  the  SIR.  This  established 
an  injector  velocity  range  for  the  model  corresponding  to  the  operating 
range  of  the  SIR. 

Thermal  similarity  is  of  no  consequence  in  this  study.  It  was 
concluded  that  with  an  injector  velocity  range  selected  as  in  this 
similarity  discussion,  the  velocity  field  in  the  model  would  relate 
well  to  that  existing  in  the  SIR.  Hence,  it  seems  to  be  worthwhile 
to  conduct  an  analog  experiment  with  water  injected  into  water. 

2 . 2  Analysis  and  Computer  Code 

A  special  case  of  the  SIR,  where  the  inflow  and  outflow  arc  on 
the  same  end  of  a  dead  ended  cylinder  was  considered  in  this  analysis, 
figure  1  depicts  the  geometry.  The  analysis  was  carried  out  using  a 
modified  form  of  TI'ACII. 

TIACII-T  is  a  form  of  the  more  general  program  TTAHI!  and  is  used 
specifically  for  computation  of  flows  in  ducts  wiMi  sudden  enlargements. 
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Modifications  were  made  to  TF.ACII-T  to  account  for  the  geometry  of  the 
SIR.  Refer  to  flosman,  et  al.  (13)  for  a  more  detailed  description  of 
TF.ACH-T. 

In  TIZACH-T  the  geometry  is  as  that  of  Figure  2,  which  shows  the 
region  of  flow  bounded  by  the  duct  walls,  the  axis  of  plane  of  symmetry, 
inlet  and  outlet  boundaries.  The  inlet  boundary  is  located  in  the  plane 
of  the  step,  where  the  conditions  of  the  incoming  fluid  are  presumed 
to  be  known;  while  the  outlet  boundary  is  positioned  sufficiently 
enough  downstream  that  the  assumption  of  fully  developed  flow  is  fairly 
accurate. 

The  inflow  is  presumed  to  be  steady  and  turbulent,  with  a  mean 

velocity  U.  .  The  diameters  of  the  large  and  small  ducts  are  designated 
in 

by  P  and  <1  respectively. 

The  predictions  of  the  problem  require  the  solution  of  the  equa¬ 
tions  of  motion  for  axial  velocity  II,  radial  velocity  V,  pressure  !’, 
together  with  the  two  equation  model  for  turbulent  closure  involving 
turbulent  kinetic  energy,  k,  and  dissipation,  a. 

The  numerical  technique  and  constants  were  retained  in  their 
entirety,  slight  changes  were  made  in  the  grid  sice.  The  only  major 
modifications  made  to  deal  with  the  model  of  Figure  1  were  changes  in 
the  boundary  conditions.  The  following  guidelines  were  used  when 
making,  changes  in  the  boundary  conditions: 

1.  Profiles  for  II,  V,  P,  k  and  i.  are  supplied  at  the  inlet  plane. 

2.  At  the  boundary  walls,  the  calculation  of  the  velocity 
component  parallel  to  the  wall,  turbulent  kinetic  energy, 
and  dissipation  are  made  using,  a  logarithmic  law  of  the  wall, 
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as  in  the  original  program. 

3.  At  the  axis  of  symmetry,  the  condition  of  vanishing 
normal  flux  is  applied  to  all  variables. 

4.  At  the  outlet  nlane,  the  normal  velocities  are  deduced 
from  the  overall  mass  continuity  requirements  and 
specification  of  a  fully  developed  velocity  profile. 

Basic  changes  to  be  made  in  the  boundary  conditions  were  two¬ 
fold: 

1.  Introduction  of  an  end  wall  at  the  desired  location 
using  the  corresponding  boundary  conditions. 

2.  Removal  of  the  boundary  condition  at  the  existing  side 
wall  and  replacement  with  an  outflow  condition. 

This  called  for  changes  in  two  parts  of  the  program  TT.ACI1-T.  Modific 
tions  were  made  in  the  main  program  to  define  the  initial  variable 
distributions  to  change  geometry  (example,  ratio  of  d  to  HI  and  to 
define  new  variables  and  changes  in  format  statemen's.  Changes  were 
also  required  in  SUBROUTlNIi  PR0M0P,  responsible  for  the  modification 
of  the  finite  difference  equation  terms  and  coefficients  at  the 
boundaries.  Changes  in  the  individual  sections  of  this  subroutine 
are  described  as  follows: 

1.  b.N'TRV  MODI):  responsible  for  modifications  for  calculation 
of  the  axial  velocity.  The  conditions  at  the  top  wall  and 
axis  of  symmetry  remain  the  same.  The  outlet  boundary  is 
replaced  with  an  end  wall;  as  shown  in  Appendix  A.  the  value 
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of  the  coefficient  of  the  finite  difference  equation  is 
calculated  and  introduced  in  this  section.  The  exist  ini’ 
side  wall  conditions  are  removed  and  replaced  by  an  outflow 
condition.  Two  outflow  conditions  were  tried.  In  one,  it 
was  assumed  that  the  exit  velocity  profile  was  fully  developed; 
in  the  other,  an  exit  velocity  profile  was  assumed  at  the 
outset . 

2.  ENTRY  MOPY:  responsible  for  modifications  for  calculation 

of  the  radial  velocity.  Side  wall  conditions  were  introduced 
in  accordance  with  the  law  of  the  wall  as  shown  in  Appendix  A. 

a.  ENTRY  MODTE:  responsible  for  modification  for  calculation 
of  turbulent  kinetic  energy.  Side  wall  conditions  were 
introduced  in  accordance  with  the  law  of  the  wall  as  in 
ENTRY  MODY. 

4.  ENTRY  MOPED:  responsible  for  modification  for  calculation 
of  dissipation  rate.  Side  wall  conditions  were  introduced 
in  accordance  with  the  law  of  the  wall  as  in  ENTRY  MOPY. 


All  other  parameters  remained  unchanged.  A  list  of  the  modified  state¬ 
ments  is  included  in  Appendix  R. 
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CHAPTER  111 

EXPERIMENTAL  PROCEDURE 

3. 1  Experimental  Method 

Qualitative  information  on  the  flow  structure  in  the  model  was 
obtained  using  a  particle  visualization  technique  similar  to  the  one 
presented  by  Ekchian  and  Hoult  (24).  Their  method  utilized  neutral 
density  particles  to  study  the  flow  structure  in  the  water  analog  of 
an  internal  combustion  engine.  Allen  and  Yerman  (25)  found  that 
neutral  density  particles  do  not  alter  the  flow  field;  quantitative 
velocity  predictions  made  from  their  technique  were  within  five  per¬ 
cent  of  pitot  tube  data.  In  this  study  alumina  particles  were  used; 
which  for  small  enough  particle  sizes  can  be  expected  to  follow  t he 
flow  as  well  as  neutral  density  particles  (26,27). 

A  two  dimensional  flow  field  was  assumed,  and  a  plane  containing 
the  axis  of  the  model  was  assumed  to  be  representative  of  the  flow- 
structure  in  the  entire  flow  field.  A  plane  so  described  was  illuminated 
at  a  high  intensity.  The  width  of  the  illuminated  zone  was  made  small 
in  comparison  with  other  dimensions  and  it  was  assumed  to  be  a  plane. 

The  particles  (assumed  to  move  along  streamlines  (26))  were  photo¬ 
graphically  recorded  as  they  passed  through  the  illuminated  plane, 
yielding  a  picture  of  the  instantaneous  flow  structure. 

3 .2  Experimental  Components 

Figure  5  is  a  schematic  diagram  of  the  water  loop.  The  major 
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component  of  the  apparatus  was  the  horizontal  cylindrical  test  chamber 
(model).  The  fluid  jet  was  injected  centrally  from  one  end  of  the 
chamber,  and  fluid  was  continuously  discharged  through  small  holes  on 
both  ends.  The  other  components  included  a  measuring  tank,  receiving 
tank,  dye  bypass,  centrifugal  pump,  valves  and  a  differential  manometer 
The  working  fluid,  water,  was  pumped  through  a  strainer  and  an 
iniector  into  the  model.  A  mercury  filled  differential  manometer  used 
across  the  iniector  nozzle  was  calibrated  to  determine  the  iniector 
velocity.  When  the  line  to  the  dye  chamber  was  opened,  water  mixed 
with  green  dye  (Merriam  0-29.70 )  and  was  injected  into  the  clear 
bath  gave  a  picture  of  the  spreading  jet. 

Outflow  from  the  model  could  lie  controlled  using  the  valves 
on  either  end  of  the  model.  Volumetric  flow  rate  was  measured  using 
a  stop  watch  (least  count  of  a  hundredth  of  a  second)  and  the  measuring 
tank.  When  dye  studies  were  performed,  the  colored  water  was  drained 
out  of  the  measuring  tank,  the  water  loop  was  open,  lor  the  particle 
visualization  studies,  i.e.  water  in  the  receiving  tank  was  seeded 
with  alumina  particles  and  the  tank  was  kept  Siirred:  water  from  the 
measuring,  tank  was  not  drained,  but  let  into  the  receiving  tank  instead 
Air  in  the  model  was  removed  by  opening,  the  air  purge  valve  on 
the  manometer  line.  When  this  valve  was  opened,  the  manometer  was 
isolated  from  the  rest  of  the  system  by  shutting  the  manometer  control 
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figure  4  is  a  picture  of  the  water  loop  as  set  up.  A  more 
detailed  description  of  the  individual  components  follows. 


3.2.1  Model  of  SIR 

The  model  was  used  in  two  basic  configurations;  the  first  one  so 
as  to  enable  dye  injection  with  a  syringe  at  various  wall  locations 
downstream  of  the  injector  plane,  the  second  one  for  particle  visuali 


The  model  consisted  of  a  visualization  cylinder  cut  out  of  a 
plexiglass  tube  (Cadco  Acrylic).  The  tube  had  an  outer  diameter  of 
4.25  inches  and  an  internal  diameter  of  3.5  inches.  A  6  inch  long 
cylindrical  tube  was  used  for  the  dye  studies  and  13.25  inch  long 
one  for  particle  vi sual  i  cat  ion,  giving  length  to  diameter  ratios  of 
1.71  ami  3 . 50  rospoc t i ve 1 y . 

Figure  5  shows  a  visualization  cylinder  as  used  for  the  dye 
visualization  studies.  A  pressure  tap  was  located  0.5  inches  from 
the  end  on  the  top  of  the  cylinder  wall,  and  numerous  other  taps, 
fitted  for  syringe  injection  wore  located  along  the  cylinder  wall. 
The  injection  fittings  were  made  by  adapting  standard  fitlings  as 
shown  in  Figure  (>.  The  glands  ensured  that  there  was  no  leakage  out 
of  the  model,  the  needle  guide  ensured  that  the  needle  (which  was 
h  inches  long,)  did  not  buckle  when  pushed  through  the  gland.  both 
ends  of  the  visualization  cvlinder  were  tapped  axially  for  twelve 
equally  spaced  10-52  screws.  Onlv  the  pressure  tan  existed  on  the 
cvlinder  used  for  the  particle  visualization  studies. 
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The  visualization  cylinder  was  used  in  conjunct  ion  with  end  caps, 
perforated  plates,  rubber  gaskets,  visualization  box  end  plates  and 
a  nozzle,  assembled  as  shown  in  Figure  7  (picture).  Figure  8  shows 
the  schematic  of  the  same. 

The  visualization  box  is  used  in  the  case  of  particle  visualization 
studies,  when  filled  with  water  it  removes  most  of  the  distortion 
in  photographs  caused  by  the  curvature  of  the  cylinder. 

Figure  0  shows  a  perforated  end  plate,  cut  0.25  inches  thick  from 
a  4.5  inch  diameter  plexiglass  rod  (Cadco  Acrylic).  The  perforations 
arc  0.125  inch  in  diameter,  set  out  in  a  0.5  inch  x  0.5  inch  matrix 
as  shown.  Aluminum  rivets,  either  blank  or  with  through  holes  fit  into 
these  perforations,  allowing  for  control  of  the  outward  mass  flux 
distribution.  The  diameters  of  the  holes  in  the  rivets  were  0.010  inch, 
0.025  inch,  0.051  inch  and  0.055  inch. 

The  end  cans  cut  from  a  4.5  inch  diameter  plexiglass  rod  were 
machined  to  shape.  They  provided  outlets  for  outflow  to  t lie  measuring 
tank. 

The  nozzle  was  machined  down  from  a  1  inch  brass  bar  stock;  it  had 
a  0.5  inch  outer  diameter  and  a  0.0095  inch  hole  diameter.  It  was 
5  inches  long  and  was  held  in  place  using  a  teflon  compression  fitting 
screwed  into  the  end  can. 

The  entire  assembly  was  held  together  by  two  sets  of  10-52  screws, 


one  set  on  each  end. 
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3.2.2  Dye  Bypass 

Figure  3  shows  a  dye  bypass  loop.  It  was  used  when  dye  injection 
through  the  nozzle  was  desired.  It  was  possible  to  divert  all  the  flow 
or  a  part  of  the  flow  through  this  bypass  where  it  picked  up  the  green 
dye  on  its  way  to  the  injector.  The  bypass  had  three  control  valves 
as  shown  in  Figure  3  and  a  dye  container  which  had  to  be  recharged 
after  each  run. 

3.2.3  Differential  Manometer 

A  Meriam  Instruments  60  inch  mercury  differential  manometer  was 
calibrated  to  give  the  injector  velocity  as  a  function  of  the  pressure 
drop  across  the  nozzle  measured  in  inches  of  mercury. 

3.2.4  Measuring  and  Receiving  Tanks 

A  plastic  measuring  tank  with  dimensions  12  inches  x  12  inches  x  12 
inches  was  used  with  a  stop  watch  to  measure  the  volumetric  flow  rate  in 
the  water  loop.  It  was  possible  to  measure  separately  the  water  flow 
rate  coming  out  of  each  end  of  the  model.  This  tank  emptied  into  the 
receiving,  tank  for  the  particle  visualization  runs,  for  dye  studies 
this  tank  was  cut  off  from  the  receiving  tank  and  all  the  water  was 
drained  out.  The  mean  water  temperature  was  found  to  be  70°F. 

The  receiving  tank  was  also  made  of  plastic,  slightly  larger  in 
size,  20  inches  x  24  inches  x  24  inches  and  normally  hel  '  110  litres  of 
water.  Alumina  particles  were  seeded  into  the  water  loop  through  this 
tank  which  was  kept  well  stirred.  A  wire  mesh  at  the  tank  outlet  acted 
as  a  strainer,  but  let  the  small  alumina  particles  through. 
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3.2.5  Forced  Circulation  Pump 

The  water  flow  through  the  loop  was  driven  by  a  pumn  located 
downstream  of  the  receiving  tank.  The  pump  manufactured  by  the 
Allis  Chalmers  Corporation  was  a  centrifugal  pump  with  a  maximum 
rpm  of  3450.  It  was  driven  by  a  0.75  IIP  close  coupled  electric  motor. 

It  could  provide  a  maximum  flow  rate  of  1.76  litres  per  second.  The 
flow  rate  through  the  system  was  controlled  by  means  of  a  throttling 
valve  on  the  pump  discharge. 

3. 2.6  Visualisation  Box 

A  rectangular  plexiglass  box  of  dimensions  5.5  inches  x  4.75  inches 
x  14  inches  was  used  to  enclose  the  visualization  cylinder  when  con¬ 
ducting  particle  visualization  studies.  The  box  was  constructed  from 
0.25  inch  thick  plexiglass  plates  glued  together  and  made  watertight 
by  using  silicon  sealant.  The  end  plates  of  this  box  fitted  into  the 
model  assembly  between  the  perforated  plate  and  the  end  plate  as  shown 
in  figure  8.  This  box  had  two  ports  on  the  top,  one  for  water  feed 
and  the  other  for  air  escape.  It  had  a  drainage  port  on  the  bottom. 

Mien  filled  with  water,  it  removed  most  of  the  optical  distortion 
caused  by  the  curvature  of  the  cylinder  and  enabled  good  photography. 
This  box  was  not  used  when  doing  dvc  visualization  studies. 

3 . 3  fighting  l.quipment 

The  light  source  consisted  of  a  115V,  in  amp  General  fleet ric 
mercurv-arc-capillary  lamp,  lamp  housing,  transformer  and  tripod. 

The  lamp  emitted  high  intensity  ultraviolet  radiation.  The  lamp  and 
housing,  were  both  air  cooled,  using  air  throttled  from  the  machine 
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shop  air  supply  at  55  psi.  The  cooling  was  monitored  by  a  f.cneral 
llectric  mercury  lamp  control  unit  which  operated  on  115'.',  100  II; 
power  supply,  obtained  using  a  brushless  inductor/alternator  made 
by  the  Leach  Co. 

The  lamp  housing  also  contained  a  system  of  parallel  slits  and 
plano-convex  condensing  lenses,  which  enlarged,  col  1  imated  and  defined 
a  thin  light  beam. 

The  light  source  and  lenses  in  a  single  housing  were  mounted  on 
a  tripod  as  shown  in  the  picture  of  Figure  10.  This  allowed  for 
raising  or  lowering  and  swivel  motion  of  the  light  housing. 

The  slits  transformed  the  light  into  a  plane  containing  the  axis 
of  the  model.  The  light  was  shielded  so  as  to  prevent  light  leaving 
from  anywhere  else  except  the  slit. 

5 .  -1  Photographic  laju i  pinont 

A  55  mm  Canon  AF.-l  camera  was  used  to  photograph  the  illuminated 
particles  through  the  transparent  visualization  box.  The  camera  was 
used  in  conjunction  with  a  12  mm  Vivitar  extension  tube  so  as  to  allow 
und i stort cd  close-up  photography.  The  camera  was  mounted  10.1  incites 
vertically  above  the  model  axis  and  the  focus  was  set  at  infinity  to 
give  a  sharp  and  clear  picture.  Provisions  were  made  for  lateral 
movement  of  the  camera,  so  the  entire  flow  field  could  be  photographed 
using  exposures  at  the  different  locations  along  the  axis  of  the  model. 
Figure  11  shows  the  camera  position  relative  to  the  light  source  and 


v  i  suii  1  i  zat  ion  box. 


si 


Anr  1  1,  1980 
id 


X  .  5  V  i  sua  1  i  :nt  ion  ‘’articles 
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s»nk  in.’  velocity.  I  he  s  i  nk  i  n  •  .  >  1  i  %  d :  rect  1  v  proportional  to 

the  stniare  of  t  he  particle  d :  ar.eter ,  md  to  t  lie  difference  bet  .eer. 
the  particle  and  fluid  density;  it  is  inversely  proportional  to  the 
fluid  viscosity  (do,.  1  hough  the  density  match  between  alumina  and 

water  is  not  too  good ,  it  is  compensated  for  by  using  very  small 
particle  sizes.  The  residence  time  of  a  fluid  element  in  the  model 
is  so  small  that  sinking  is  almost  non-existent.  Thus,  such  particles 
arc  expected  to  map  the  flow  well. 

The  traces  of  the  particles  as  seen  on  the  photographs  were  pro¬ 
jections  of  t  lie*  i  r  paths  into  a  plane  containing  the  axis  of  the  model. 

The  particles  were  seeded  into  the  water  loop  through  the  receivin'',  tank 
which  was  kept  stirred. 

X.(>  Photography 

Optimum  photography  was  obtained  using  a  Kodak  Tri  X  film  (TX - 1  da -.10 1 


which  was  black  and  white  film  exposed  and  developed  at  100  ASA.  The 
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optimum  camera  setting  was  obtained  at  a  shutter  speed  of  1/15  second 
and  an  aperture  f8.0  when  the  light  slit  was  8  to  10  mm  wide. 

The  negatives  were  used  to  project  an  image  on  a  screen  using 
a  film  strip  projector. 


i 
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CILAPTllK  IV 

PRESENTATION  AND  I  NTEUPRETAT  I  ON  OP'  RP.SUI.TS 

4.1  Dye  Visualization  Experiments 

Two  techniques  were  used  for  dye  visualization  studies  in  the 
model.  In  one  case,  a  dye-water  mixture  was  injected  into  a  clear 
bath  in  the  model,  through  the  nozzle.  In  the  second  case,  a  syringe 
was  used  to  inject  pure  dye  into  a  clear  hath  in  the  model,  through 
glands  fitted  to  the  model  wall. 

4.1.1  Nozzle  Pvc  Injection 

Test  runs  were  made  using  inlet  velocities  ranging  from  3  m/s 
to  IS  m/s.  These  velocities  corresponded  to  those  in  the  SIR 
applying  the  similarity  conditions  of  Chapter  II.  The  corresponding 
inlet  Reynolds  number  range  was  from  4170  to  405.30.  According  to 
McNaughton  and  Sinclair  (IS),  for  inlet  Reynold  numbers  greater  than 
3000  for  1 iquid- l iquid  jets  in  short  cylindrical  vessels,  fully  turbulent 
icts  are  observed.  Indeed,  these  dye  visualization  studies  did  reveal 
a  turbulent  jet.  comprised  of  a  forward  moving  conical  jet  region  and 
a  slower  moving  backflow  region  which  surrounded  the  jet.  No  laminar 
length  was  observed,  and  the  jet  mushroomed  at  the  point  of  injection. 

The  dye  used  was  neutral  density  and  was  assumed  to  follow  the 
streamlines.  The  technique  however  could  not  provide  more  detailed 
information  in  this  case  due  to  the  diffusivity  associated  with 
turbulent  flows.  High  flow  velocities  caused  the  entire  bath  to  cloud 
in  a  fraction  of  a  second,  further  limiting  the  effectiveness  of  the 
techniuue.  Ihese  experiments  established  that  the  jet  was  turbulent 
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and  that  recirculation  was  predominant. 

•1.1.2  Syringe  Dye  Injection 

This  technique  was  intended  to  help  in  gaining  insight  into  the 
reattachment  lengths  of  the  turbulent  iet.  It  was  also  hoped  that  this 
method  would  provide  a  means  of  isolating  particular  regions  in  the 
flow  that  were  of  interest.  A  similar  technique  was  used  by  Back 
and  Roschke  (16),  in  which  dye  was  bled  into  the  shear  layer  between 
the  jet  and  reverse  flow-  region,  through  holes  in  a  circular  channel 
wall  unstream  of  a  sudden  expansion.  Their  studies  were  limited  to  a 
Reynold  number  range  from  20  to  4200. 

I>yc  was  injected  into  the  bath  at  various  locations  downstream 
of  the  noetic.  The  syringe  was  inserted  to  different  depths  and  the 
dye  released.  By  releasing  dye  at  a  given  axial  location,  but  varying 
radial  locations,  the  shear  layer  was  located  at  the  point  above  which 
the  dye  moved  in  the  reverse  direction  and  below  which  the  dye  moved 
with  the  iet  flow.  When  dye  was  injected  into  the  shear  layer,  it  was 
expected  to  follow  it  until  it  reattached  at  the*  cylinder  wall  at  some 
downstream  location.  Unfortunately,  due  to  the  high  turbulence 
intensity  and  fluctuating  flow  pattern  such  observation  was  not  clear. 
At  lower  Reynolds  numbers  (about  52 (if*  to  -00n)  it  was  possible  to 
infer  that  the  flow  did  not  reattach  to  the  cylinder  walls:  when  it 
approached  the  end  of  the  cylinder  the  dye  seemed  to  turn  and  a  part 
of  it  would  flow  back  along  the  cylinder  wall.  It  must  be  pointed  out 
that  due  to  the  clouding  of  the  dye  at  the  end  of  the  cylinder  and  due 
to  the  fluctuating  location  of  the  shear  layer,  such  inferences  arc  not 
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very  obvious  and  could  be  questioned.  These  inferences  however  seem 
to  be  in  agreement  with  Back  and  Roschke's  predictions  that  reattachment 
occurs  about  9  to  10  step  heights  downstream  of  injection  for  turbulent 
jets;  in  this  case  the  cylinder  was  only  3.5  step  heights  long. 

This  technique  established  the  fluctuating  nature  of  the  flow  and 
that  numerous  eddies  existed  in  the  reverse  flow  regions.  These  eddies 
were  isolated  since  they  would  contain  dye  injected  into  them  for  a 
small  length  of  time  until  fluctuations  caused  the  eddy  to  disappear. 

This  time  period  was  not  long  enough  for  detailed  observation.  The 
need  was  felt  for  obtaining  an  overall  picture  of  the  flow  structure 
and  to  observe  how  this  overall  structure  fluctuated. 

4,2  Particle  Visualization  F.xper  iment  s 

Illuminated  alumina  particles  were  used  to  visualize  the  flow- 
field  in  the  model  with  a  length  to  diameter  ratio  of  3.5.  The 
experiment  was  run  in  the  jet  Reynolds  number  range  similar  to  that 
m  the  dye  experiments.  Outflow  at  both  ends  of  the  model  was  varied. 

Figure  12  shows  the  typical  flow  pattern  observed.  The  flow- 
consisted  of  a  spreading  jet,  a  strong  recirculating  back  flow  near 
the  wall,  and  numerous  eddies  which  .appeared  and  disappeared  in  less 
than  1/10  of  a  second.  These  eddies  caused  the  jet  flow  to  bend 
around  them,  their  appearing  and  disappearing  causing  a  side  to  side 
whipping  motion  of  the  jet.  It  seemed  that  the  jet  whipped  from  one 
wall  to  another  in  an  almost  regular  manner.  This  whipping  was  stronger 
for  higher  Reynolds  numbers  but  was  confined  to  a  region  near  the  extreme 
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end  (away  from  the  injector)  of  the  cylinder.  For  small  Reynolds 
numbers,  the  whipping  was  weaker  but  was  prevalent  in  the  entire  flow 
field;  in  this  case  the  jet  would  whip  in  a  serpentine  manner.  These 
observat ions  are  not  obvious  from  Figures  13,  14  and  15  as  these 
figures  represent  instantaneous  flow  pictures. 

Figures  13,  14  and  15  provide  the  contrast  between  varying  outflow 
at  either  end.  All  three  pictures  represent  an  inlet  Reynolds  number 
of  7057.  The  three  frames  in  each  picture  represent  three  segments 
of  the  model  each  photographed  separately  and  at  different  instants 
of  time. 

Figure  15  shows  the  case  where  all  the  outflow  is  on  the  extreme 
end  of  the  model.  The  jet  entrainment  near  the  injector  end  can  be 
seen.  Recirculation  is  not  very  strong  compared  to  the  case  of 
Figure  15,  in  which  outflow  is  all  on  t he  injector  end.  Alumina 
particles  tend  to  accumulate  at  the  extreme  end  where  they  arc  dis¬ 
charged.  One  can  observe  a  lot  of  eddy  structure  in  the  central  segment, 
surrounding  the  high  velocity  jet  core  which  is  wrinkled  and  broken. 

In  Figure  15,  not  many  particles  are  visible  in  the  backflow 
region  along  the  walls  due  to  the  high  velocities  existing  in  this 
case  which  makes  it  difficult  to  record  the  particle  motion  at  the  given 
level  of  illumination.  In  this  case  a  stagnation  region  is  observed 
at  the  extreme  end  and  particles  can  be  seen  to  turn  away  from  the 
direction  of  jet  flow  into  the  reverse  flow  region.  In  this  case  the 
number  of  eddies  existing  are  not  as  many  as  in  Figure  15.  let  entrain¬ 


ment  is  evident. 
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I  *!'.•  1  •>  Picture  of’  Plow  field  at  Reynolds  Number  70S? 
(outflow  on  injector  end) 


*! 
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Figure  14  shows  the  case  where  the  outflow  is  equally  divided 
between  the  two  ends  of  the  model.  It  does  seem  that  this  case  is 
between  those  of  Figures  13  and  15.  Particle  accumulation  is  absent 
at  the  extreme  end  and  the  reverse  flow  velocity  near  the  wall 
becomes  large  only  for  the  near  injector  region.  Eddies  and  wrinkled 
jet  structure  arc  predominant  in  the  central  segment. 

It  is  concluded  that  the  outflow  configuration  is  of  major 
consequence  in  determining  the  flow  structure  in  the  model. 

The  effect  of  higher  Reynolds  numbers  was  studied.  Figure  16 
shows  the  picture  of  the  flow  structure  for  a  Reynolds  number  of 
32315,  it  is  the  case  with  all  the  outflow  on  the  injector  end  and 
is  compared  with  the  similar  case  (Figure  15)  for  a  Reynolds  number  of 
7057.  The  higher  Reynolds  number  case  shows  a  highly  turbulent  jet 
with  a  high  velocity  backflow.  There  is  high  jet  entrainment  in  the 
near  noccle  region  and  a  resulting  high  spreading  rate  of  the  jet.  A 
lot  of  small  eddy  structure  is  visible  in  the  extreme  end  region. 

This  is  in  contrast  with  the  lower  Reynolds  number  cases  where  eddy 
structure  was  predominant  throughout  the  length  of  the  model.  Figure 
17  shows  the  same  outflow  configuration  for  a  jet  Reynolds  number  of 
1800.5.  The  overall  structure  remains  the  same,  small  eddy  structure  is 
ovserved  even  in  the  central  segment . 

4 . 3  Results  from  the  Modified  TFA(il  (Computer  ('ode 

TF.ACI1  was  run  for  the  model  using  a  water  analog  with  outflow 
onlv  on  the  injector  end,  the  extreme  end  was  a  region  of  stagnation. 
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i.R-  1~  Picture  of  Plow  field  at  Reynolds  Xumhcr  1S005 


(outflow  on  injector  end! 
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The  outlet  velocity  profile  was  specified  in  two  ways: 

1.  Assumed  a  specified  uniform  exit  velocity  profile. 

2.  Assumes  a  fully  developed  velocity  profile  at  the 
exit  plane. 

Figure  18  shows  that  regardless  of  how  the  exit  velocity  profile 
is  specified,  the  corresponding  axial  velocity  profiles  are  within 
+5"  of  each  other  at  an  x/L  of  just  0.065.  As  x/L  increases  the 
two  profiles  approach  each  other  more  closely.  It  is  concluded  that 
despite  such  big  differences  in  the  specified  exit  velocity  profile, 
the  downstream  velocity  profiles  are  not  very  different:  the  choice 
of  outlet  condition  is  not  too  important  as  long  as  continuity  is 
satisfied.  This  comparison  was  made  for  a  case  with  an  injector  di¬ 
ameter  one  third  the  size  of  the  cylinder  diameter,  it  is  assumed 
that  the  same  trend  continues  for  smaller  injector  diameters,  such  as 
that  in  the  model. 

In  all  subsequent  runs,  a  fully  love  loped  velocity  profile  was 
assumed  at  the  exit  plane.  This  provided  for  a  more  continuous 
variation  in  velocity  profile  from  one  axial  location  to  the  next. 

Figure  19  shows  the  comparison  between  axial  velocity  profiles  in 
the  pipe  step  case  and  the  modified  case  for  the  model. 

The  residues  from  iteration  to  iteration  show  good  and  stable 
convergence.  The  values  of  the  residues  for  tie  II  momentum,  V  momentum 
and  continuity  become  very  small,  as  does  that  for  turbulent  kinetic 
energy.  The  dissipation  residue  though  small  appears  to  be  cm  ey  ■; • 
at  a  slower  rate.  These  runs  were  made  using  the  same  under  re!  iv.it.  n 
factors  that  were  used  in  the  original  program.  4  he  value--  of  the  •Via' 
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residues  were  only  slightly  higher  than  in  the  original  test  case. 

Figure  20  shows  the  velocity  vectors  at  some  of  the  grid  points 
for  the  case  where  the  injector  Reynolds  number  25590,  comparable  to 
that  of  Figure  16.  Appendix  C  tabulates  data  For  this  case  The  figure 
shows  a  strong  iet  core  with  slight  entrainment.  The  backflow  is 
strong  near  the  wall  and  eddies  appear  to  exist,  sue  as  in  the  stagnation 
region  near  the  extreme  end.  Results  obtained  at  both  higher  and  lower 
Reynolds  numbers  show  the  same  trend.  F.ntrainment  is  very  great  in  the 
injector  vicinity  as  can  be  seen  by  the  velocity  vector  li  in  Figure  20. 
Another  interesting  trend  observed  using  the  program  at  very  high  Reynolds 
number  is  the  appearance  of  additional  small  eddies.  It  is  important 
to  realize  that  all  eddies  will  not  show  up  on  the  computer  results 
since  time  averaging  causes  some  of  the  regularly  fluctuating  eddies 
to  disappear. 
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CHAI’TLR  V 

SUMMARY  AN!)  CONCLUSIONS 


5.1  Summary 

The  objective  of  this  study  was  to  investigate  the  fluid  flow 
structure  of  a  jet  induced  mixing  process  in  a  SIR,  where  the  iet 
is  directed  along  the  cylinder  axis. 

A  water  model  was  constructed  which  satisfied  the  conditions  of 
similarity.  The  analogy  leads  to  the  conclusions  that  the  flow  pattern 
observed  in  the  model  is  representative  of  that  in  the  SIR  and  that  the 
velocities  in  the  model  are  related  to  those  in  the  SIR. 

The  water  model  was  a  horizontal  cylindrical  vessel,  containing 
a  water  bath  at  room  temperature.  A  water  iet  was  injected  centrally 
from  one  end  of  the  model  and  was  continuously  discharged  through  small 
holes  on  both  ends.  The  water  flowed  in  a  closed  loop  for  particle 
visualization  studies  and  in  an  open  loon  for  dye  visualization.  Runs 
were  made  at  various  inlet  Reynolds  numbers.  The  onflow  on  both  ends 
was  also  varied. 

Two  techniques  were  used  for  the  dye  visualization  studies  in 
the  model.  Nozzle  dye  injection  was  employed  to  study  the  spreading 
behavior  and  the  turbulent  nature  of  the  flow  field.  Syringe  dye 
injection  through  glands  in  the  model  wall  provided  for  a  study  of 
iet  rcattachment  and  eddy  structure. 

The  particle  visualization  entailed  photographing  the  generated 
flow  patterns  which  were  made  vis ib le  by  illuminating  small  alumina 
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particles  placed  in  the  water  loop.  The  particles  were  visually  ob¬ 
served  in  the  model  when  they  passed  through  a  narrow  zone  of  high 
intensity  illumination.  The  narrow  zone  of  light  approximated  a 
plane,  containing  the  cylinder  axis  and  point  of  injiction. 

A  modified  form  of  the  TEACH  computer  code  was  used  to  make 
analytic  predictions  of  velocity ,  pressure,  turbulent  kinetic  energy, 
and  dissipation  fields  in  the  specific  case  of  a  dead  ended  model 
with  all  outflow  in  the  plane  of  injection. 

5.2  Results  and  Conclusions 

For  a  jet  induced  mixing  process  in  a  SIR,  the  major  observations 
and  conclusions  of  this  investigation  are: 

1.  The  flow  consisted  of  a  spreading  turbulent  iet  with 
a  strong  region  of  reverse  flow  near  the  walls.  Eddy 
structure,  predominant  in  the  reverse  flow  region, 
indicated  high  turbulence  intensity. 

2.  Fluctuating  eddies  caused  the  jet  core  flow  to  bend 
around  them,  producing  a  side  to  side  whipping  motion 
of  the  jet.  The  flow  structure  was  highly  fluctuating. 

5.  The  jet  did  not  reattach  to  the  cylinder  walls. 

4.  At  high  injector  Reynolds  numbers,  tiie  whipping  was 
strong  but  confined  to  the  extreme  end  of  the  model, 
where  eddy  structure  was  concentrated.  The  reverse 
flow  along  the  wall  was  strong.  .Jet  entrainment  near 
the  injector  was  high  causing  a  rapid  spreading. 
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For  lower  Reynolds  numbers,  the  whipping  motion  was 
weak.  Fddies  were  prevalent  throughout  the  flow  field 
causing  a  serpentine  whipping.  The  reverse  flow  was 
weaker  and  the  entrainment  (hence  spreading)  less. 

Outflow  configuration  was  felt  to  have  a  major  effect  on 
the  flow  structure.  The  effect  of  higher  outflow  on  the 
injector  end  was  to  produce  higher  reverse  flow  veloci¬ 
ties  near  the  walls,  eddy  structure  became  less  prevalent, 
and  a  region  of  stagnation  began  to  form  at  the  extreme 
end.  Increased  outflow  at  the  extreme  end  caused  smaller 
reverse  flow  velocities,  eddy  structure  became  predominant 
especially  in  the  central  segment  of  the  cylinder  length, 
the  stagnation  region  at  the  extreme  end  disappeared. 

Results  obtained  using  the  modified  form  of  TFACIi  arc  time 
averaged  flow  field  variables,  for  a  case  with  no  outflow 
on  the  extreme  end.  Results  showed  a  strong  turbulent  jet 
with  a  high  velocity  reverse  flow.  A  stagnation  region  was 
encountered  at  the  extreme  end,  where  the  streamlines  of  the 
jet  bend  into  the  region  of  reverse  flow.  Fddy  structure 
was  indicated. 

Comparison  between  visual  and  analytic  trends  was  not 
feasible.  This  was  due  to  the  fact  that  a  time  averaged 


visual  structure  could  not  be  constructed. 
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5.3  Recommendations  for  Further  Study 

It  would  be  of  great  value  to  record  the  flow  patterns  on  movies 
or  video  tape.  This  would  provide  for  a  better  understanding  of  the 
fluctuating  structure  and  provide  a  base  for  time  averaged  visual 
observations. 

Velocity  measurements  should  be  made  to  permit  a  thorough 
evaluation  of  the  analytic  model.  The  laser  doppler  velocimeter 
would  seem  to  be  the  proper  choice  of  instrument  for  such  turbulent 
recirculating  flow  configurations. 

The  analytic  model  has  potential  for  further  modifications.  It 
could  be  extended  to  account  for  outflow  at  both  ends. 

It  is  also  suggested  that  a  wider  range  of  Reynolds  number  and 
length  to  diameter  ratio  he  studied. 


53 


April  1,  1980 
SL/CHW : lei 


REFERENCES 


1.  Thomson,  F.  X.,  "Investigation  of  a  Jet-Induced  Mixing  Heat 
Transfer  Process,"  M.S.  Thesis,  The  Pennsylvania  State  University, 
1974. 

2.  Becker,  11.  A.,  Hottel,  H.  C.  and  Williams,  0.  C. ,  "Mixing  and 
Flow  in  Ducted  Turbulent  Jets,"  Ninth  Symposium  (International) 
on  Combustion,  Academic  Press,  London,  1963,  pp.  7-20. 

3.  Barchion,  M.  and  Curtet ,  R. ,  "Some  Details  of  the  Structure  of 
an  Axisymmetric  Confined  Jet  with  Backflow,"  Transactions  of 

the  ASMH,  Journal  of  Basic  Engineering,  December  1964,  pp.  777-78" 

4.  Craya,  A.  and  Curtet,  R. ,  "On  the  Spreading  of  a  Confined  Jet," 

(in  French)  Comptcs-rendus  Academic  des  Sciences,  Paris,  Vol. 

241,  1955,  pp.  621-622. 


5.  Curtet,  R. ,  "Confined  .Jets  and  Recirculation  Phenomena  with  Cold 
Air,"  Combustion  and  Flame,  London,  Vol.  2,  No.  4,  December  1958. 

6.  Curtet,  R. ,  "On  the  Flow  of  an  Enclosed  .Jet,”  Publications 
Scicntif'iqucs  et  Techniques  du  Ministcre  de  l’Air,  Paris, 

No.  359,  1  60. 

7.  Abramovich,  0.  N.  ,  The  Theory  of  Turbulent  Jets,  The  M.I.T.  Press, 

1963.  . 

8.  Sehlichting,  II.,  Boundary  Layer  Theory,  McC.raw-lli  1 1  Co.,  196S. 

9.  Hill,  P.  C.,  "Turbulent  Jets  in  Ducted  Streams,"  Journal  of  Fluid 
Mechanics,  Vol.  22,  Part  1,  1965,  pp.  161-1S6. 

10.  Korney,  P.  J.,  E'aeth,  0.  M.  and  Olson,  D.  It.,  "Penetration  Charac¬ 
teristics  of  a  Submerged  Steam  Jet,"  AIChE  Journal,  Vol.  IS,  No.  3, 
1972,  pp.  548-554. 

11.  Weimcr,  J.  C.  ,  laeth,  0  M.  and  Olson,  D.  R.  ,  "Penetration  of 
Vapor  Jets  Submerged  in  Subcooled  Liquids,”  AIChE  Journal,  Vol .  19, 
No.  3,  1973,  pp.  552-557. 


Avery,  J.  F.  and  Faeth,  0.  M.  ,  "Combustion  of  a  Submerged  Caseous 
Oxidizer  Jet  in  a  Liquid  Metal,"  Fifteenth  Symposium  (International) 
on  Combustion,  Academic  Press,  New  York,  1974,  pp.  501-512. 


54 


April  1,  1980 
SL/CHW:lcl 


13.  Gosman,  A.  D. ,  Khalil,  E.  E.  and  Whitelaw,  J.  H. ,  "The  Calculation 
of  Two-Dimensional  Turbulent  Recirculating  Flows,"  First  Inter¬ 
national  Symposium  on  Turbulent  Shear  Flows,  Springer-Ver lat 
Berlin  Heidelberg,  1979,  pp.  237-255. 

14.  Gosman,  A.  D. ,  Khalil,  E.  E.  and  Launder,  B.  E.,  "Computation  of 
Flow  in  Ducts  with  Sudden  Enlargements,"  Turbulent  Recirculating 
Flows  Workshop,  The  Pennsylvania  State  University,  July  1975. 

15.  Launder,  B.  E.  and  Spalding,  D.  B.  ,  Mathematical  Models  of  Turbulence, 
Academic  Press,  London,  1972. 

16.  Back,  L.  H.  and  Roschke,  E.  J.,  "Shear-Layer  Flow  Regimes  and  Wave 
Instabilities  and  Reattachment  Lengths  Downstream  of  an  Abrupt 
Circular  Channel  Expansion,"  Journal  of  Applied  Mechanics,  Vol.  94, 
1972,  pp.  677-681. 

17.  Reynolds,  A.  J. ,  "Observations  of  a  liquid-liquid  jet,"  Journal  of 
Fluid  Mechanics,  Vol.  14,  1962,  pp.  552-556. 

18.  McNaughton ,  K.  H.  and  Sinclair,  C.  G.  ,  "Submerged  Jets  in  Short 
Cylindrical  Flow  Vessels,"  Journal  of  Fluid  Mechanics,  Vol.  25, 

1966,  pp.  367-375. 

19.  Crow,  S.  C.  and  Champagne,  F.  II.,  "Orderly  Structure  in  Jet  Turbu¬ 
lence,"  Journal  of  Fluid  Mechanics,  Vol.  48,  1971,  pp.  547-591. 

20.  White,  W.  N. ,  Darooka,  D.  K.  and  Sogin,  11.  N. ,  "On  Recirculating 
Flow  in  a  Dead-End  Vessel,"  Transactions  of  the  ASME,  Journal  of 
Fluids  Engineering,  March  1975,  pp.  119-120. 

21.  Willis,  D.  A.  and  Meyer,  W.  E.,  "A  Water  Model  Study  of  Swirling 
Flow  in  a  Stratified  Charge  Engine,"  Department  Report,  Mechanical 
Engineering,  The  Pennsylvania  State  University,  1965. 

22.  Bird,  R.  B.  ,  Stewart,  W.  E.  and  I.ightfoot,  E.  N.,  Transport  Phenomena, 
John  Wiley  and  Sons  Inc.,  2nd  edition,  1962. 

23.  Prandtl,  L.  and  Tictjcn,  0.  J.,  Applied  Hydro-  and  Aero-dynamics, 

Dover  Publications  Inc.,  1957. 

24.  Ekchian,  A.  and  lloult,  I).  P.  ,  "Flow  Visualization  Study  of  the 
Intake  Process  of  an  Internal  Combustion  Engine,"  SAE  Technical 
Paper  Scries,  Detroit,  February  1979. 

25.  Allen,  M.  and  Yerman,  A.  J.,  "Neutral  Density  Beads  for  Flow 
Visualization,"  Section  4,  Symposium  of  Flow  Visualization, 
Presentation  Summaries,  A SMI.  Annual  Meeting,  Chapter  4,  I960. 


55 


April  1,  1980 
SL/CHW: lcl 


Merzkirch,  W. ,  Flow  Visualization,  Academic  Press,  New  York,  1074. 

Wood,  C.  J.,  "Visualization  of  an  Incompressible  Wake  with  Base 
Bleed,”  .Journal  of  Fluid  Mechanics,  Vol.  29,  1967,  pp.  259-272. 

Fulmer,  R.  D.  and  Wirtz,  D.  P. ,  "Measurement  of  Individual  Particle 
Velocities  in  Simulated  Rocket  Exhaust,"  AIAA  Journal,  Vol.  5, 

No.  8,  1965,  pp.  1506-1508. 


56 


April  1,  1980 
SL/CHW: lcl 


APPENDIX  A 

CALCULATION  OF  FINITE  0 1  FF FRENCH  EQUATION 
COEFFICIENTS  FOR  NEW  SILL  WALL 

Refer  to  Figure  21  for  the  nodal  configuration  at  the  introduced 
extreme  end  wall.  The  boomerangs  show  had  data  is  stored. 

Consider  the  node  (IS.d)  where  J  can  vary  from  2  to  21,  the  flux 
term  (general)  is  defined  as; 


i  =  i  'u  i  (ft 
'  e  c  c 


r 

1  ;A  3x  e 


For  the  axial  velocity  momentum  equation, 

r.  3ll  I 


1  =  cu  U 

e  e  e 


u  3xe 
e 


since  u  is  verv  small  and  negligible  in  the  near  wall  region 

O'  ‘ 

i  =0,  therefore  Au  =  0. 

•  e  e 

S ini  lari  lv  for  the  \'  momentum; 


1  =  i'll  v 

•  e  e  e 


v 

amYsT 


?v 


r  g_‘_! 
v  3.x '  e 


V  ...viiro 

wa  1  1 


e 

Y(15) 


v  am  ir.) 


therefore 


\v  =  0,  and  S'  =  S  -  -v 
e  p  p  Ax 

these  changes  in  the  source  term  S  are  made  in  accordance  with  the  law 
of  the  wall,  as  used  in  the  original  program.  Similar  treatment  is  used 


for  $  -  k ,  •  . 
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APPENDIX  B 

MODIFICATIONS  i N  Till:  COMP1HT.R  CODE 


The  following  changes  were  made  in  the  computer  code: 

Changes  in  the  MAIN  PROGRAM 
CHAPTER  I 

new  values  were  assigned  in  the  following  statements,  as  shown 
N.J  =  22 

R LARCH  =  0.0445 
ALTOL  =  7.0  *  R LARGE 

Statements  numbered  45  and  44  were  removed 
The  following  statements  were  added: 

RAT  =  RSMALL/RLARGE 

UOUT  =  -UIN*RAT**2/ (1 . 0- RAT* *21 

TF.OIJT  =  TURBIN*  UOUT* *2 

HPOUT  =  TEOUT*  *1.5/  ( A  LAM  DA  *  R  LA  RG  E ) 

CHAPTER  II 

The  following  statements  were  included: 

DO  201  J  =  JSTPI,  N.1M1 
TE(l,.n  =  TEOUT 
U(  2,  .11  =  UOUT 
201  ED(1  ,.I)  =  EDOirr 

DO  205  I  =  2,  N I M I 
205  yPLUSN ( I )  =  11.0 

Statement  number  117  was  replaced  by; 

DO  657  .1  =  2,  .1STEP 


Changes  in  ‘‘UBROUTl NE  PROMOD 


ENTRY  MODI! 

C - TOP  WALL 

No  changes  required. 

C - SIDE  WALL 

DO  215  .1  =  2.N.IM1 
215  AEINIMI,.!)  =  0.0 

C - SYMMETRY  AXIS 

No  changes  required. 

C - OUTLET 

ARD!  NT=0. 0 
El 0W=0. 0 

DO  204  ,T=. JSTPI  ,  N.IM1 

ARDEN  =  0 , 5*  (.DEN  ( 2  ,.!  1 +DLN  ( 5  ,.1)  1  *R  ( .11  *SNS  (.1) 
ARDENT  =  ARDENT  ♦  ARDEN 
204  FLOW = FLOW  +  ARDEN*ll(5  ,.I1 

II  INC  =  !  E  LOW  I N  *  FLOW) /ARDENT 
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no  205  .J-.JSTP1  ,N.IM1 
05  11(2, J)  -  l)(5,.J)  -UINC 

Rr.TURN 

LNTRY  MOOV 

c - sim:  wall 

COTLRM  =  CMU**0. 25 
XP  =XU (N I ) -X (NIM1 1 
I  =  N I M 1 

no  510  .1=2  ,N.1M1 

SQRTK=SORT(P .  5*  (TL (I  ,.J)  +TL  ( I , J-l ) ) ) 
nrs'v=  o.5*(ni:N(i  ,j)+nr.N(i  ,.j-id 
X P LUSA=  0 . 5 *  ( XPLUSL  (-7 )  +.XP l.( ISE  (.) - 1 )  ) 

IP  (XPI.USA.LL.il. 631  00  TO  SI  1 

TMULT-!)i:\V*(:n'n.RM*.SQRTK*CAPPA/AL0n(nL0n*XPHISA) 
no  TO  312 

11  i'MWI.T-YISCO.S/XP 

12  !  Mil  M  I  =-TMULT*Y(  I  ,.n 
10  AL  (  I  ,.J  I  =0.  0 

TAUL.f  11-TAUL(21 
TAUL  (N.J )  -TAIIL  l  .X.T'tl  I 

C - TOP  WALL 

Nn  changes  required. 

C - -  SYMMI  TRY  AXIS 

No  change  required. 

RLTURN 

LNTRY  MODI’ 

RLTURN 

LNTRY  MOOT 
RLTURN 

LNTRY  MOIVIL 

C - TOP  WALL 

No  changes  required 

C - SI  PL.  WALL 

XP=X IK.NI )  - X  (NIM1 ) 

I  =N  I M 1 

no  (,20  .1=2,  N. JM 1 
ni.NY=l)LN(I  ,.J) 

SQRTK-SQRTrn:  ( 1  ,.!)) 

Y0L=R  (.1 1  +SNS  (.1 1  *SLW ( 1 1 

XPI.iISi;(J)  =  nLNV*SQRTK*a)TLR'rXP/\'I.S<:OS 

0LN(T)U=0 . 58  (ABS  (I  AIM.  (.1+  1 1 '  Y  ( I  ,.I*  1  ))♦  MIS  (TAUL  (.11  *Y  ( I  )  ) /XP 
nvnx=  (  ( V  (  I  ,.n+Y(  I  ,.T+1  1  *Y(  I  4  1  ,.!)  *  V  ( I  ♦  I  ,.J+  1  11/  I .  -  (Y(  I  ,.n*V(i  ,.]+n 
♦  VI  1-1  ,.1  )  »Y(  !-  1  ,.M  i  )/-l.O)/si  W  (  I  ) 

ULNI’I  S~ULN  1  f  ,.M  -VISC  I  ,.1 1  '*U\  UX**^ 

OLNl  I  ,.!)=ULNRI  S  +  ULNCOII 
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1 1:  ( Xl’l.USF.  ( J ) .  LL:.  .11.631  HO  TO  621 

lHTFRM=nFN(l  ,J )  *  (CMl)**0.75)  *SQRTK*ALOn(F.LOr,*xrLUSE(.J)  )  /  (CAPP.\*XP1 
no  TO  622 

621  CONTINUE 

!>ITFR'I=DFN  (I ,J) *  (CM1I**0. 75)  * sqrtk*xplusf  (J )  /XV 

622  CONTI  NUT. 

su(i,.i)=su(i  ,.j+sum>(i  ,.n+r.r.N(i,j)*voL 
SP  ( I  ,J)  =  SP(I,J)+SPKD(I  ,.I)-UITKRM*VOL 
620  AP.(I,.n=0.0 

c - SVMMTTRV  AXIS 

No  changes  required 
RF.TURN 

FNTRY  MODFP 

C - TOP  IV A  L L 

No  changes  required 

c . '  sim;  kali. 

Xr=Xl)(NI)-X(NlMl) 

I=N1M1 

TFRM= (CMU**0. 75) ? (CAPPA'XP) 

NMM2=N.J  2 
DO  720  <J=2 , N.JM2 

SU(  I  ,.J )  =ORFAT*TFRM*TF  ( I  ,J)**1.5 
720  SI'fl  ,.n  =  -C, RF. AT 

C - SYMMFTRY  AXIS 

No  changes  required. 

RFTUR.N 

FNI) 
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APPENDIX  C 

AXIAL  AND  RADIAL  VELOCITIES  AS  CALCULATED  BY  THE  MODIFIED  TEACH 
PROGRAM  CORRESPONDING  TO  A  REYNOLDS  NUMBER  OF  2S590. 
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55E+00 

0. 

4 

5. 

2Gi; 

+  00 

5 , 

251 ; 

+  00 

5. 

2oi; 

♦  00 

5. 

29! 

+  00 

5, 

24 1; 

+  00 

5. 

,  66F.+  00 

4  , 

25F.+00 

0 . 1 

5 

9. 

491 

+  00 

9, 

251: 

+  00 

8. 

9ir: 

+  00 

8. 

54 1; 

+  00 

8. 

.  15E  +  00 

8. 

.  10E+00 

4. 

94E+00 

0.' 

9. 

901. 

+00 

9. 

ss  r; 

+  00 

9. 

851: 

+  00 

9. 

77i; 

+  00 

9. 

.  651. 

+  00 

9. 

.571':+ 00 

5, 

59F.+00 

0. 

X  = 

0, 

,  0G9 

0. 

.  089 

0. 

,  112 

0. 

,  159 

0. 

,  172 

0. 

.212 

0, 

.260  0 

.51 
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tabu:  11 

RADIAL  V  L  LOG  I T I  OS  AT  GRID  POINTS  AS  CALCULATED  BY  MODIFIED  TEACH 


I 

T 

3 

4 

5 

6 

7 

8 

->  A 

0. 

0 

0. 

0 

0. 

0 

0 

.0 

0. 

0 

0. 

0 

0. 

.0 

21 

1. 

36  E 

-03 

1. 

19E 

03 

7. 

17E 

-04 

4 

15F. 

-05 

-7. 

1SE 

-04 

-1  . 

5  OF. 

-03 

_  n 

.2  OF: 

■03 

20 

3 . 

25E 

-03 

3 . 

12E 

03 

04F. 

-03 

6 

.  39E 

-04 

-8. 

97  F. 

-04 

_  -> 

45F. 

-03 

-3. 

.  84  E 

-03 

19 

5. 

26  E 

-03 

5. 

23E 

03 

3. 

49  E 

-03 

1 

3SF. 

-03 

-9. 

02E 

-04 

-3. 

19E 

-03 

-5. 

03 

IS 

- 

16E 

-03 

*7 
/  , 

23E 

03 

4. 

84  E 

-03 

2 

,  08E 

-03 

-8. 

89E 

-04 

-  3 . 

87E 

-03 

-6. 

.  49E 

03 

1? 

8. 

"9E 

-03 

8. 

91E 

03 

5. 

93E 

-03 

"> 

.611- 

-03 

-9, 

8  OF. 

-04 

-4. 

59E 

-03 

-  7 . 

.  75E 

03 

16 

9. 

9'E 

-03 

1 . 

OOE 

02 

6. 

581: 

-03 

.  84E 

-03 

-1 . 

291; 

-03 

-5. 

47E 

-03 

-9 . 

,  05E 

03 

15 

1. 

06  E 

-02 

1. 

04E 

02 

6. 

69  E 

-03 

2 

,  67E 

-03 

-1 . 

92E 

-05 

-6. 

59E 

-03 

-1 . 

,  05E 

02 

14 

1 . 

041: 

-02 

9. 

S4E 

03 

6. 

20E 

-03 

2 

.  02E 

-05 

_2 

99E 

-03 

-8. 

OOE 

-03 

-1  . 

.  23E 

02 

13 

<1 , 

2SE 

-03 

8. 

5 1 E 

03 

5. 

21E 

-03 

8 

,  391: 

-04 

-4. 

65  E 

-03 

-1  . 

OIF. 

-02 

-1  . 

.  4SE 

02 

12 

/  , 

461: 

-03 

6. 

951. 

03 

4. 

08E 

-03 

-8 

.851: 

-04 

04F. 

-03 

-1  . 

29  E 

-02 

-1 

.  73E 

02 

11 

5. 

311 

-03 

6. 

18E 

03 

3 . 

24  E 

-03 

-  3 

.131: 

-03 

-1 . 

04  E 

-02 

-1  . 

66E 

-02 

_  1 

.  OSE 

02 

10 

3 . 

13E 

-03 

6. 

23E 

03 

-> 

19E 

-03 

-6 

,  1 3  E 

-03 

-1 . 

48E 

-02 

_  -> 

14E 

-02 

_  -1 

.  521: 

02 

9 

-4, 

25  E 

-04 

-1. 

42E 

03 

-4. 

65  E 

-03 

-1 

.  16E 

-02 

_  -> 

05F. 

-02 

__  -> 

72E 

-02 

-  3 , 

.  OOE 

02 

8 

-1 . 

761; 

-03 

-6. 

131: 

03 

-1 . 

281; 

-02 

_  -> 

.  OOE 

-02 

_  -> 

671: 

-02 

-3 . 

111; 

-02 

-3 

.  1  OF. 

-02 

■? 

-1. 

77  E 

-03 

-6. 

561: 

-03 

-1 , 

331: 

-02 

_  "A 

.Oil: 

-02 

-2. 

4  81: 

-02 

_  -> 

66F. 

-02 

_  -> 

.  55E 

-02 

6 

2 

OOE 

-04 

-6, 

561; 

-05 

-  I . 

74  E 

-03 

-4 

.  87E 

-03 

-  7 , 

93  E 

-03 

-9. 

09  E 

-03 

-7, 

.  94  E 

-03 

5 

8. 

45E 

-04 

3. 

47E 

-03 

7 , 

29E 

-03 

1 

.  201; 

-02 

1, 

61 E 

-02 

1  . 

S2E 

-02 

1  , 

.  80E 

-02 

4 

1 . 

47E 

-03 

1 . 

72E 

-03 

54  E 

-03 

4 

.  8 1 E 

-03 

8 

95  F. 

-03 

1 . 

45E 

-02 

1  , 

.  97E 

-02 

3 

3 , 

35  E 

-05 

1 , 

59  E 

-04 

1  , 

HE 

-04 

2 

.  4 1 E 

-05 

-3 , 

781: 

-05 

-1  . 

08  E 

-05 

1  . 

.  96E 

-04 

T 

0. 

0 

0. 

0 

0, 

0 

0 

.0 

0. 

0 

0. 

0 

0. 

.0 

X  = 

0. 

.  003 

0. 

.  009 

0 

.016 

0, 

.024 

0. 

.  034 

0. 

046 

0 

.061 
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TABLE  II  (cont.) 


9 

10 

11 

12 

0.0 

0.0 

0.0 

0.0 

-2. 77E-03 

-3.25E-03 

-4. 10E-03 

-5.34E-03 

-4 . 94E-03 

-5 . 84F.-03 

-7.46E-03 

-9.67E-03 

-6.28E-03 

-8.07E-03 

- 1 . 04E-02 

-1.33E-02 

-8.51E-03 

-1.00E-02 

-1.30E-02 

-1.61E-02 

-1.01E-02 

-1. 18E-02 

-1.52E-02 

-1.81E-02 

-1.  17E-02 

- 1 . 34E-02 

- 1 . 73E-02 

- 1 . 9 1 E  -  0  2 

-1.34E-02 

- 1 . 49E-02 

-1.92E-02 

-1.92E-02 

-1.52E-02 

-1.63E-02 

-2. 13E-02 

-1 .S2E-02 

-1.74E-02 

-1 . 78E-02 

-2.35E-02 

-1.59E-02 

-2. 00E-02 

-1.94E-02 

-2.60E-02 

-1.25E-02 

-2. 32E-02 

-2.90E-02 

-2.90E-02 

-7.20E-03 

-2. 70E-02 

-2.23E-02 

-3.20E-02 

-3. 901: -04 

-3.06E-02 

-2.29E-02 

-3.28E-02 

7.151-:- 02 

-3.00E-02 

-2. 14E-02 

-2.68E-02 

1. 17E-02 

-2. 13E-02 

-1.44E-02 

-1.49E-02 

1.25E-02 

-4.97E-03 

- 1 . 42E-03 

-7.78E-04 

1.44E-02 

1.68E-02 

1.52E-02 

1.35E-02 

1.86E-02 

2.32E-02 

2. 4 IE- 02 

2.35E-02 

2.34E-02 

7. 16E-04 

1.58E-03 

3.09E-03 

4.77E-03 

0.0 

0.0 

0.0 

0.0 

0.078 

0.099 

0.124 

0.  154 

13 

14 

15 

16 

0. 

0 

0. 

0 

0 

0 

0. 

0 

-9. 

211:- 

03 

-1  . 

93E- 

02 

5 

1 4  E  - 

02 

0. 

0 

-1. 

75E- 

02 

-3 . 

89E- 

02 

1 

02E- 

01 

0. 

0 

-2. 

50E- 

02 

-5. 

88E- 

02 

1 

52F.- 

01 

0. 

0 

-3. 

14E- 

02 

-r 

83E- 

02 

2 

00E- 

01 

0. 

0 

-3 . 

66E- 

02 

-9 . 

65E- 

02 

2 

46F.- 

01 

0. 

0 

-4. 

03E- 

•02 

-1  . 

12E- 

■01 

2 

89E- 

-01 

0. 

0 

-4. 

25E- 

•02 

-1. 

25E- 

-01 

3. 

28E- 

-01 

0. 

0 

-4. 

55E- 

-02 

-1  . 

33E- 

-01 

3 . 

62E- 

-01 

0. 

0 

-4. 

27E- 

-02 

-1 . 

55E- 

-01 

3. 

90E  • 

-01 

0. 

0 

-4. 

,  10E- 

-02 

-1. 

,  29E- 

-01 

4. 

.  12E- 

-01 

0. 

0 

-  3 . 

.811:- 

-02 

-1 . 

,131;- 

-01 

4, 

27E- 

-01 

0. 

,0 

-3 , 

.  4SE- 

-02 

-8 

.  47  E 

-02 

4 

34  E 

-01 

0 

.0 

-3 

.  1 4  E  ■ 

-02 

-4 

.  261: 

-02 

4 

321:  ■ 

-01 

0 

.0 

-3 

.  19E- 

-02 

1 

.4  71: 

-02 

4 

201: 

-01 

0 

.0 

-3 

.  77E 

-02 

S 

.  031: 

-02 

5 

961: 

-01 

0 

.0 

-  3 

.  43E 

-02 

1 

.  33E 

-01 

3 

551: 

-01 

0 

.0 

-1 

.  81E 

-02 

1 

.611: 

-01 

2 

961: 

-01 

0 

.0 

-> 

.  031: 

-05 

1 

.591: 

-01 

2 

17E 

-01 

0 

.  0 

1 

.  85E 

-03 

9 

.  75E 

-02 

1 

15E 

-01 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.0 

0 

.190 

0 

.234 

0 

.  286 

0 

.  33 
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